N87- 22214 


ESTIMATION OF SQUEEZE-FILM DAMPING AND INERTIAL COEFFICIENTS 

FROM EXPERIMENTAL FREE-DECAY DATA 

J.B. Roberts 
University of Sussex 
Brighton, England 

R. Holmes 

University of Southampton 
Southampton, England 

P.J. Mason 

Chelsea College of Aeronautical and Automobile Engineering 
Shoreham, Sussex, England 


This paper describes the results obtained from an experimental programme 
concerned with a parametric identification of the damping and inertial coefficients 
of a cylindrical squeeze-film bearing, through an analysis of transient response 
data. The results enable the operating range for which a linear model of the 
squeeze -film is appropriate to be determined. Comparisons are made between the 
estimated coefficients and theoretical predictions. Presentation is by courtesy of 
the Council of the Institution of Mechanical Engineers, London. 

INTRODUCTION 

Accompanying the development of modern machinery has been an increasing demand 
for higher running speeds. As a result critical speeds have been encountered before 
the desired running speed of the machine has been reached. The consequences of this 
can result in failure of associated components or, in some cases, an inability to 
reach the desired operating speed. The present work is concerned with the 

squeeze-film damper, which is proving very successful in mitigating these problems. 
A manageable and realistic model for a squeeze-film damper can be derived by 
applying linearisation techniques to the oil -film forces, which are obtained by 
solving the Reynolds equation. This leads to a representation of the dynamic 
behaviour in terms of damping coefficients. The ability to provide damping is a 
feature of this device but there is no capacity to provide linear stiffness as the 
latter depends on journal rotation. 

The conventional representation of a squeeze -film in terms of damping coeffi- 
cients has the attraction that it is very simple to incorporate these coefficients 
into a discrete mathematical model of a rotor-bearing system. This approach 
implicitly assumes that inertial forces within the oil -film are negligible. 
However, theoretical investigations by a number of workers [1-11] suggest that 
inertial forces can be very significant in squeeze -films. Indeed, this is evident if 
one considers the "gap Reynolds number” 
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where p is the density of the lubricant, o> is the frequency of vibration, c is the 
radial clearance in the squeeze -film damper and rj is the viscosity of the lubricant. 
In many practical applications R e is of order one, or greater. For such values one 
cannot expect conventional lubrication theory, based on the Reynolds equation (ie. 
on the assumption that R e « 1), to give an accurate representation. 

In a linearised approach, inertial forces can be incorporated into the model 
through the introduction of "acceleration coefficients” in addition to the normal 
damping coefficients. However, experimental estimates of acceleration coefficients 
(or "hydrodynamic mass" effects), in geometries typical of many squeeze -films, do 
not appear in the literature; indeed, inertial effects have not been taken into 
account at all in many previous comparisons between theoretical and experimental 
dynamic behaviour. This may explain, at least partly, why it has proved so diffi- 
cult to obtain satisfactory agreement between theoretical and experimental values of 
the damping coefficients. 

In the present investigation a transient testing technique has been developed 
and used to obtain the dynamic characteristics of a "short” squeeze-film bearing, 
with a geometry typical of that currently adopted in engineering applications and 
with Reynolds numbers in the range 0.5 < R e < 1.5. The experimental results are 
processed using a parametric identification technique to yield estimates of the 
dynamic coefficients and these are compared with predictions from existing theory. 
Limitations of the theory are highlighted which indicate promising avenues for 
further research. 
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NOTATION 

direct fluid damping coefficients, for the r and s directions, 
respectively. 

structural damping coefficients, for the r and s directions, 
respectively. 

non-dimensional, direct fluid film damping coefficients, for the r 
and s directions, respectively. 

radial clearance between the journal and the bearing. 

direct fluid film inertial coefficients, for the r and s direc- 
tions , respectively. 

initial displacement 

non-dimensional initial displacement (see equation (27)). 
time domain response function (see equation 17). 
shaft stiffnesses, in the r and s directions, respectively, 
land length of the bearing, 
effective, first mode mechanical mass. 

hydrodynamic mass (c rr in the radial direction, c ss in the 
transverse direction). 
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hydrodynamic film force. 

journal displacement in the radial direction (i.e. in the 
direction of the attitude line). 

radius of the journal. 

Reynolds number (=pwc 2 /h) 

journal displacement in the transverse direction (i.e. in a 
direction perpendicular to the attitude line). 

Lime. 

displacement . 

non-dimensional displacement, 
measured , f ree - decay record . 

frequency response function, defined by equation (23). 
frequency response function, defined by equation (23). 
damping parameter, defined in equations (3). 
value of y, in the absence of a fluid film, 
non -dimensional ised y (see equations (10)). 
stiffness parameter, defined by equations (3). 
value of 8, in the absence of a fluid film, 
non -dimensional ised 8 (see equations (10). 

critical damping factors, in the case of no fluid film, in the r 
and s direction, respectively. 

eccentricity ratio. 

static eccentricity ratio. 

absolute viscosity of the squeeze-film fluid, 
density of the squeeze-film fluid, 
sampling interval, 
frequency of vibration. 

natural frequencies of undamped vibration, in the radial and 
transverse directions, respectively. 

non-dimensional time (= o> 0 t). 
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THEORETICAL TREATMENT 


The equations of motion 


The system consists of a mass of finite size (the squeeze-film journal) 
attached to the centre of a beam, which is built-in to rigid supports at each end 
and which plays the role of a conventional retainer spring. If the central mass is 
given some initial .transverse displacement, achieved by a force applied at the 
centre of the journal, and then released, the first mode of vibration will be 
dominant in the subsequent motion. Thus to a close approximation, the motion can be 
described in terms of two second order equations of motion. 


In the present investigation, the squeeze-film bearing was run under full-film 
conditions, that is with no cavitation. In these circumstances, it can be shown, 
theoretically, that no cross -damping terms appear. 

We thus write, for radial r and transverse s, displacements 
(m + c rr )r + (b° r + b rr )r + k r r = O 

(m + c ss )s + (b° s + b ss )s + k s s = O ... (1) 

Here c denotes inertial coefficients and b damping coefficients in the 
squeeze -film, b° denotes structural damping and k structural stiffness emanating 
from the spring beam. m is the effective, first mode mechanical mass. Both 
equations can be written in the standard form 


x + yx + Sx = 0, 

where y and 6 are constants. For the radial direction: 


( 2 ) 


_ b° r + b rr 
in + Cpp 

... O) 


in 4- Cp* p 

and similarly for the transverse direction. 


From a free decay test, one can determine y and 6, by using a parametric 
identification technique. Suppose that a free decay test is carried out in the 
radial direction, in the absence of a fluid film in the squeeze -film bearing. Then 
b rr = c rr = O, and we obtain the coefficients 

b° r 

Vo - ~ 2 £r w o 

k • - * ( ^ ) 

6 ° ~ in” = Wr ° 2 

where u) r0 is the natural frequency of undamped vibration and £ r is the damping 
factor. On combining equations (3) and (4) one obtains 


b r r — m 


r 6 0 

y — - y 0 


Crr — m ( " 1 ) 


... (5) 
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for the damping and inertial coefficients, where 


m 


®o 


m can be determined from measurements of k r and from a free decay test without 
fluid in the squeeze-film bearing (giving y 0 and 6 0 ). By combining this information 
with estimates of y and 6, from free decay with fluid in the squeeze-film bearing, 
one can estimate b rr and c rr , by using equation (5). Similarly, measurements of 
free decay in the transverse direction can be used to yield estimates of b S9 and 
c 88 . 


In practice the structural damping is very small and one finds u ro and u 80 by 
singly measuring the frequency of vibration, without fluid in the bearing. 

Before processing the decay curves it is convenient to non-dimensionalise the 
equation of motion (equation (2)). If u 0 is the frequency of undamped vibration, 
without fluid, then a convenient non-dimensional time is 

... (7) 

can be non -dimensional ised by dividing 

Thus, 

... ( 8 ) 


0 ... (9) 


. . . ( 10 ) 


are non-dimensional coefficients and differentiation is now with respect to r. On 
substituting these equations into equation (5) we obtain 


t = w 0 t = 6 0 “t 

Also, the displacement x(t), 
by the initial displacement, x(0). 


x (t) = 

x( O ) 

Equation ( 2 ) can then be recast as 
X + y* X + 6* X = 

where 


y* =- y 
u>o 


_ y 


6 a 

6* = - , = ^ 

Wq 2 6 C 


brr = mwro O*/®* - Yo*] 
Crr = m(^ - 1) 

and similarly for b 33 and c S3 . 


The damping coefficients 

The conventional approach to evaluating the coefficients b rr and b ss , for a 
full squeeze-film is to use the Reynolds equation as the basis of the calculation. 
For the general case it is necessary to solve this equation numerically, but simple, 
asymptotic results can be derived for 

(a) the short bearing: fi/R -» o 

(b) the long bearing: £/R -» oo. 


297 


The short bearing theory gives the following results [12] 
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( 12 ) 


The corresponding long -bearing results are as follows[13] 
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( 13 ) 


As pointed out earlier, the use of the Reynolds equation implicitly assumes 
that inertial forces within the fluid film are negligible. Theoretical studies of 
the influence of fluid inertia on the damping coefficients of a squeeze-film bearing 
have been undertaken by Tichy[3] and San Andres and Vance[7]. Their results 
indicate that, for the particular geometry of bearing studied in this investigation, 
and for the range of frequencies of oscillation studied, the influence of fluid 
inertia on the damping coefficients is negligible. 


The inertial coefficients 


A linearised approach to the evaluation of squeeze -film fluid forces allows a 
separate evaluation of the inertial coefficients, which arise from the effect of 
journal acceleration. 


Smith[l] has shown that for a very short full -film bearing ( ««R) these 
coefficients are given by, 


Cp r 


ffpRfi 3 2 1 

‘ a - «./? 


- 1]} 


c ss 


irpR£3 
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[1 - (1 - e, 


> 2 )^]} 


... ( 14 ) 


and c rs = c sr = 0. In the special case of zero static eccentricity (e 0 
limiting operation performed on equation (14) gives 


0), a 


c r r = c ss rrpRg 3 (e 0 = o) 


... (15) 


12 C 


This agrees with the result obtained by Fritz for very short bearings, in the 
concentric case [2]. 
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For a long, full -film bearing ( fi»R) the appropriate expressions for the 
coefficients are as follows [l]s 


ffpR3g , 2 - k 

c rr - c ss — { g [1 - (1 - e 0 2 )*} 


(16) 


and c rs = c 3r =0. In the special case where e 0 = 0 (concentric operation) equation 
(16) reduces to 


rrpR3fi . 

c r r - c S3 - - — (e 0 = 0) 


... (17) 


which is a result first derived by Stokes [14]. 

It is interesting to note that, for e 0 = o, the ratio of c rr for the long bearing 
case to c rr for the short bearing case, is from equations (15) and (17). 


The damping coefficients from the long and short bearing theories, for the 
concentric case, are in exactly the same ratio (see equations (12) and (13)). It 
follows that, for short bearings (fi « R), the damping and inertial coefficients 
are, according to the short bearing theory, considerably less than those predicted 
from the long bearing theory. 

Recently Szeri et al [9] have presented, graphically, numerical values for the 
inertial coefficients, c rr and c ss , for squeeze - fi lms with finite values of fi/R, in 
the range 0.2 < £/R < 4.0. They found it necessary to introduce an approximation, 
based on the assumption that fi/R is small, in their analysis but have indicated that 
their results should be more accurate than a full short -bearing approximation, 
provided that i/R Is small. Their numerical results are in virtually exact agree- 
ment with Smith's short -bearing results (equation (14)) for g/ R <0.5. 

2.4 Parametric identification 

It was shown in section 2.1 that the damping and inertial coefficients can be 
related to the non-dimensional parameters y* and 6*, which occur in the second-order 
linear equation of motion given by equation (9). From the experiment to be des- 
cribed one can obtain a free -decay curve - i.e. X(t) versus time. The problem is 
then to find, the values of y* and 6* for which the solution to equation (9) gives a 
"best fit" to the experimental observations. 

This is a problem in parametric identification, on which considerable litera- 
ture exists (for example see Ref. [15]). Of the various available techniques, we 
have here selected the sequential method of Detchmady and Stidhar[16], since this 
enables estimates to be obtained from a knowledge of a single measured state 
variable (here the displacement of the journal versus time - either r( t ) or s(t). 
Details of the algorithm are given in the Appendix. 

The algorithm operates on a discretely sampled record of the decay curve. 
Suppose the experimental values are of y(iAr) ( i=0,l,2, . . . ), where at is the 
sampling interval and r is given by equation (7). The data is conveniently scaled so 
that the start displacement y(0), is unity. By sequentially processing the data, 
the algorithm generates least -square estimates of the state vectors X(t) and x(t). 
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and also least-square estimates of y* and 6*, at the times iAr. These estimates 
"track" the experimental data and should converge as time increases - i.e. the 
estimated state vector X(t) should approach the measured response, y(t), and the 
estimated values of y* and 8* should approach constant values. 

To obtain improved estimates of y* and 6*, the algorithm can be applied in an 
iterative manner. In the first iteration, values of y* and 8* are guessed and used 
to start the sequential estimation computation . The algorithm will give estimates 
of y* and 8*, at the end of the data sequence; i.e. at time t* , where r* is the 
time of the last data sample. These estimates can now be used, in place of the 
initial guesses, as a start to the second iteration. By repeating the iteration a 
number of times, the estimates of y* and 6*, at time T m , should converge to constant 
values . 


Prior to using this procedure on real decay data it was tested thoroughly on 
simulated data, from which it could be concluded that the algorithm was an efficient 
and useful method, for the present application. 

Memory effects 

The use of damping and inertial coefficients is based on the assumption that 
the fluid film forces depend only on the instantaneous velocity and acceleration of 
the journal. Although the coefficient approach li simple to incorporate into a 
discrete mathematical model of a rotor-bearing system, there are two serious 
objections which can be raised, concerning its validity: 

(i) Implicit in the method is the assumption that the velocity and accelera- 
tion are linearly independent variables, so far as the fluid film is 
concerned. This poses conceptual difficulties - e.g. how can the accel- 
eration be varied whilst the velocity is held constant? 

(ii) No allowance is made for "memory" effects, which can be expected when the 

bearing is running under cavitation conditions. Even in the case of a 

completely non-cavitated bearing, considered in the experimental work 
reported here, memory effects may be significant, for sufficiently high 
frequency motion, due to the visco-elastic properties of typical lubri- 
cants (e.g. see Ref. [3]). 


Considering, for example, the case of radial motion only, a general linear form 
for the relationship between the hydrodynamic force P r , and the motion, r(t), is as 
follows: 



h(t - r) r(T) dr 


(18) 


Where h( ) is a time domain impulse response function. 


The use of equation (18) allows a generalisation of the coefficient representa- 
tion discussed earlier. To demonstrate this, consider the simple case of harmonic 
motion 


r(T) = Ae iUT 


(19) 


On substituting this motion into equation ( 18 ) one obtains 
P r (t) = Ae iwt [cx(u)) + i/ 3 (u>)] 


( 20 ) 
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where 



cos uu du 


sin uu du 


... ( 21 ) 


Equation (20) may be compared with the corresponding result obtained by the co- 
efficient approach, which states that, for a squeeze -film bearing 

P r (t) = b rr r + c rr r 

Combining equations (19) and (22) one has 

P r (t) = Alwt ( -o)2 c rr + iw b rr ] 

Equations (20) and (23) are identical if 

<x(u) as -c rr u 2 

/3(<» 0 = b rr w 

However, the integral representation of equation (18) allows an arbitrary 
variation of the "coefficients" with frequency, whereas, according to the coeffi- 
cient approach, the damping and inertial coefficients are necessarily independent of 
frequency. 

It. remains to be tested by experiment whether, over a frequency range of 
practical concern, the frequency independent coefficient approach gives a satisfac- 
tory approximation, or whether there is a significant "memory effect", with the 
result that the coefficients must be treated as frequency dependent parameters. In 
the latter case, an integral representation, such as that given by equation (18) is 
more appropriate than the coefficient representation. 

DESIGN OP EXPERIMENT 

An outline drawing of the general arrangement of the rig is shown in Fig. l and 
a photograph of the rig is shown in Fig. 2. A non -rotating journal is contained 
within the circular bearing and is supported by a beam of circular cross-section 
which provides a stiffness in parallel with the squeeze-film. A cross -noct.ional 
view of the journal and beam assembly is shown in Fig. 3. The journal is heat 
shrunk along its entire contact length with the beam, and the beam is rigidly 
clamped at both ends. In the experimental work, three different beams were used, of 
varying stiffness. 

Adjustment of the static equilibrium position in the horizontal direction was 
achieved by providing a machined channel in which the bearing housing could slide. 
Care was taken to ensure that the journal was accurately aligned with respect to the 
bearing. The alignment could be adjusted by moving the position of the beam end 
supports, using shims. Angular misalignment could be effectively eliminated by 
ensuring that the distance through which the journal could be moved, within the 
bearing, was maximised. 


. . . ( 22 ) 
. . . (23) 


. .. (24) 
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The bearing consisted of two plain lands separated by a central circumferential 
groove. Lubricant was supplied from a pump, through top and bottom feed holes and 
distributed around the bearing by the groove. No end seals were fitted and the 
lubricant was free to discharge into a reservoir prior to recirculation. By 
applying an adequate supply pressure to the inlet oil, full lubricant film condi- 
tions were maintained with no cavitation. 

The experimental technique consisted of pulling back the journal, across the 
clearance circle, to a known position by a length of wire looped over the core of a 
solenoid. The journal was released by actuating the solenoid. Capacitive probes 
then transmitted the transient decay to a microprocessor based data -acquis it ion 
system. This allowed the dynamic characteristics (i.e. the mass, stiffness and 
damping of the squeeze- film, support beam and journal mass) to be evaluated and 
compared with simple linear theory. 

EXPERIMENTAL RESULTS 

In all the tests reported here the line of centres of journal and bearing, ie 
the radial direction, was horizontal when the journal was in its static equilibrium 
position. In the radial tests the journal was pulled out radially a further initial 
displacement, d, and released. Processing of the results from these tests enabled 
estimates of the damping coefficient, b rr , and the inertial coefficient c rr , to be 
derived. In the transverse tests, the journal was given an initial transverse 
displacement, d, (perpendicular to the line of centres) and released. From these 
tests, estimates of the damping coefficient, b 3S , and the inertial coefficient, c ss , 
could be derived. 

It is convenient, henceforth, to refer to a non-dimensional initial displace- 
ment, d*, defined by 

_ initial displacement (d) 

- radial clearance (c) ' ' 
Tests without fluid in the bearing 

For each of the three beams available, tests were carried out, in both the 
radial and transverse directions, with no fluid in the squeeze-film bearing. Here 
the damping is very small, and is structural in origin; it follows that the measured 
natural frequency of oscillation is, in these circumstances, a very close approxima- 
tion to the undamped natural frequency 

By applying the parametric identification procedure to the results, estimates 
of the undamped natural frequencies, w 0 , and the structural damping factor, £, were 
obtained. For a given shaft these values were found to differ slightly, in the 
radial and transverse directions (the maximum difference was about 5%), and results 
from each direction were averaged. Table 1 summarises the results obtained from 
these tests. 


Tests with fluid in the bearing 

A series of tests was carried out, with Tellus RIO as a lubricant in the 
squeeze-film bearing, and with a sufficient supply pressure to ensure that full-film 
conditions were maintained throughout (i.e. no cavitation). Measured decay curves, 
in both radial and transverse directions, were obtained for static eccentricity 
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ratios of 0, 0.1, 0.2, 0.3, 0.4 and 0.5 and for initial displacements, d*, of 0.2 
and 0.4. The decay curves were normalised, in every case, to give an initial 
displacement of unity. 

If the journal -bearing system behaves linearly, then Ute normalised decay 
curves should be independent of the initial displacement, d*, if other parameters 
are Kept constant . Thus by performing a series of decay tests, with differing 
initial displacements, and comparing normalised decay curves, one can assess the 
range within the bearing clearance circle, for which linear conditions prevail. 

At e 0 = 0.0 and c 0 = 0.3, the system behaved linearly (to a close approxi- 
mation), in both directions, for start amplitudes up to d* = 0.4. At the highest 
static eccentricity ratio studied, e 0 = 0.5, a good collapse of the normalised decay 
curves was still obtained in the transverse direction for d* = 0.2 and 0.4, whereas, 
in the radial direction there was distinct evidence of non-linearity, for d* = 0.4. 
It can be concluded that there is a fairly wide range of journal displacement 
position, within the clearance circle, for which a linear mathematical representa- 
tion is reasonable. 

Figs. 4 <a) to (f) show a set of experimental, normalised decay curves, 
obtained with the squeeze-film journal mounted on beam l . Hk; results cover the 
static eccentricity range e D = 0.0 to 0.5, and relate to both radial and transverse 
tests. Similar series of results were obtained for beams 2 and 3. In general, for 
each beam, e 0 value, and chosen direction, results were obtained for d* = 0.2 and 
0.4; where these collapsed reasonably well they were averaged to produce curves 
such as those shown in Fig. 5. Where non-linearity was indicated by a lack of 
collapse ( generally at e 0 = 0.4 and e 0 = 0.5, in the radial direction), the result 
for d* =0.2 only was used. 

A comparison between the results for e Q = 0.0 for beam 1, in the radial and 
transverse directions, (see Figs. 5(a) and (b)) indicates some degree of asymmetry 
in the journal bearing configuration. For exanqple, the second, positive overshoot in 
the radial direction is appreciably less than that observed in the transverse 
direction. The reason for this asymmetry is not clear, but may be due to the 
geometry of the oil -feed arrangement (fluid was supplied at the top and bottom of 
the central circumferential groove). The corresponding results for beams 2 and 3 
indicated that the asymmetry was much less marked at higher frequencies of oscilla- 
tion. 

Beam stiffness results 

To enable estimates of the damping and inertial coefficients to be derived from 
the free decay data, it is necessary to know the value of the effective, first mode 
mass, m. This value will be related to the actual mass of the journal, together 
with the mass of the beam, and so will vary, depending upon which beam is used in 
the experiments. 

In the face of various uncertainties regarding the precise end conditions of 
the beams, it war; decided to evaluate m, for each beam, from a knowledge of the 
natural frequency of oscillation, w 0 . and the beam stiffness. Thus 



will give an estimate of m, if k is the effective, first mode stiffness. 
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The beam stiffness, k, was determined experimentally, by pulling the 
journal, radially, with a known static force, and measuring the resulting radial 
journal displacement. Table 1 shows the results obtained for the effective masses 
of the journal and beams . 

COMPARISON BETWEEN THEORY AND EXPERIMENT 
Parameter estimation 

As a first stage in the analysis of the decay data, the parameter estimation 
procedure described in section 2 was applied to each decay curve. This yielded 
estimates of the parameters y* and 6* , in the linear, second order model given by 
equation ( 9 ) . 

Pigs. 6 show typical results of applying the estimation method to a particular 
decay curve. Here the experimental decay curve of Pig 6(b) was obtained for beam 1, 
in the radial direction, and with the journal initially concentric (e 0 = 0) (see 
also Pig. 5(a)). The iterative technique, described in section 2, was used to 
refine the estimates of y* and 6*, denoted y* and 6*, respectively. Pig. 6(a) shows 
the variation of y* and 6*, with time (measured in units of r = w 0 t) during the 
fifth iteration; at this iteration stage convergence is achieved, as evidenced by 
the fact that the final estimates in the cycle are equal to the initial estimates 
(y* =• 0.601, 6* = 0.662). Fig. 6(b) shows a corresponding comparison (for the fifth 
iteration, again) between the estimated displacement state variable, £(t), and the 
experimental decay curve; this shows that the estimated state "tracks” the experi- 
mental curve extremely well. A better idea of the degree of fit achieved can be 
obtained by comparing the experimental curve with the theoretical curve, found by 
using the final parameter estimates; this comparison is shown in Fig. 6(c). 

The excellent degree of fit obtained in Fig. 6(c) can only be obtained, of 
course, by allowing both the parameters y* and 6* to "float". If one assumes that 
the squeeze -film produces only a damping effect then it is necessary to set 6* = 1 
and to obtain a best fit by allowing only y* to vary. This can be achieved, using 
the same parametric identification procedure as before, but setting the initial, 
off-diagonal elements of the P matrix to zero (see Appendix); this has the effect 
of "locking" the 6* parameter to its initially set value, with the result that 
optimisation is sought with respect to the y* parameter alone. Fig. 6(d) shows a 
comparison between the experimental decay curve and the best -fit theoretical curve, 
with 6* = 1.0 and y* optimised (y* = 0.912). A comparison between Figs. 6(c) and 
(d) reveals that the effect of deviations of 6* from unity (due to inertial effects 
in the squeeze - fi lm ) is very significant and that a very poor fit to the data is 
achieved by assuming that only damping is present in the squeeze-film. Similar 
comparisons have been made with other decay curves and these lead to a similar 
conclusion. 


Damping and inertial coefficients 

Once y* and 6* have been estimated from a particular decay curve, then the 
damping and inertial coefficients may be found. 

To present the damping results it is convenient to introduce the non-dimen- 
sional coefficients per land. 
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Brr = brr 2^fi3 

_ q3 

Bss " bss 2fFjRp 

According to the short -bearing theory, (equation (12)), 


and 


Bpr 


(1 + 2e 0 2) 
Tl'"e 0 T) B /2 


Bss - (T 




. . . ( 26 ) 


... (27) 


In the experimental rig the following values apply: 

R = 0.06 m, c = 2.54 x 10" 4 m, fi = 0.012 m, T) = 22 x 10"3 Ns/m2 , 

Using these values in equations (26), and equations (11) one obtains an equation for 
B rr (and also B ss ) of the form 

B rr = k-j Y *-^ (i = 1,2,3) ... (28) 

where 

y*cor = y* - y 0 *$* ... (29) 

is the damping parameter, corrected for the effect of structural damping and ki is a 
non-dimensional constant, dependent on the beam used. The appropriate k^ values are 
given in Table 1 

The relationship between the present experimental values of B rr and B ss , and 
the corresponding theoretical values, according to both long and short -bearing 
theories, is shown in Figs. 7(a) and (b). It is evident that the experimental 
values lie much closer to the short -bearing theoretical curve. At zero eccentricity 
there is a factor of 300 between the two theoretical values, whereas the experimen- 
tal values are only a factor of about 1.5 higher than the short -bearing theoretical 
result. There is little indication of any "memory effect", due to changes in the 
natural frequency of oscillation. 

The inertial coefficients, c rr and c ss > defined in section 2, have the physical 
significance of hydrodynamic masses. Thus mn , the hydrodynamic mass, is given by 

m H = c rr (radially) 

= c 3S (transversely) 

Figs. 8(a) and (b) show the variations of mn with static eccentricity ratio, 
e 0 , for the radial and transverse directions, respectively. Here the experimental 
estimates of m^ are compared with both long and short bearing theoretical values. 
The significant feature here is the magnitude of the inertial effect. In both 
directions, the hydrodynamic mass is an order of magnitude greater than the short 
bearing theoretical prediction. This is rather surprising in view of the small fi/R 
ratio perlaining in the experimental rig (0.2). For this value of fl/R the results 
of Szeri et al[9] lead one to expect that the short -bearing theory should give a 
reasonably accurate estimate of mn* 
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CONCLUSIONS 


Damping coefficients for both the radial and transverse directions, agreed 
reasonably well with the short -bearing theoretical results, although the experimen- 
tal values were generally higher than the theoretical values. The variation of 
damping coefficients with static eccentricity ratio, e 0 , was very similar to the 
variation predicted by short -bearing theory. Thus, in the radial direction there 
was a marked increase in damping coefficient with e 0 , whereas in the transverse 
direction this effect was much less significant. 

The experimentally -determined inertial coefficients (or hydrodynamic masses;) 
were generally much higher than the theoretical values given by the short -bearing 
theory - typically an order of magnitude higher. 

The experimentally determined damping and inertial coefficients, for the three 
shafts, were found to collapse fairly well, when plotted against static eccentricity 
ratio. This is a strong indication that, at least over the frequency range studied 
here, "memory effects" within the squeeze film are not significant. 

APPENDIX 

Parametric Identification of Free Decay Data 
General theory 

Consider a dynamic system defined by the following differential equation of 
motion: 


* = g(x) . . . (Al) 

Here x is an n -vector containing the states, x^ , X 2 , ...,x n , of the system and 
g(x) is an~n-vector function. If the system is stable, and is released from some 
initial condition, x(O), then a transient response will result, with the motion 
decaying to zero. Suppose that observations of the output, or response, of the 
system are made during the time interval 0 < t < T. An observation vector, y(t) 
will be defined by 

y(t) = h(x) + (observation error) ... (A2) 


where y is an m-vector output and h is an m-voctor function. Here the (observa- 
tion error) term accounts for the fact that the output observation is of limited 
precision, due to quantisation errors in A/D conversion, electrical noise, etc. The 
estimation problem is to estimate the state vector x(T) from the observation 
vector y(t), measured in the interval 0<t<T-i.e. to find the vector x(T), say, 
which corresponds to a "best fit" to the observations, and is consistent with 
equation (Al). 

A "best fit" is most conveniently achieved in a least-square sense. Suppose 
that the following residual errors are defined: 

ei(t) = y(t ) - h(x*) . . . (A3) 

e 2 (t) = x* - g(x*) 
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. .. ( A4 ) 



Here x* is a "nominal” trajectory - i.e. a possible time history of x. From 
the residual errors one can form the integral 
r T 

1 = [ ll ei(t) ll + ll e2(t) ll ldt ••• < A5 > 

o J Q 2 W 2 

where j | ' | j q 2 and | | • | | w 2 are suitabl y defined quasi -norms. 

Suppose that I is minimised when x*(t ) = x(t); a least -squares estimate of x(T) is 
then x( T ) . ~ 


In practice, y(t) is usually measured at equi -spaced times, ti x iAt (i = 
0,1,2,...). It is then convenient to use a recursive algorithm, which will 
generate sequential estimates of x(t), at times t, . It has been shown by Detchmendy 
and Stidhar[16], using the method of invariant imbedding, that x(t) can be generated 
sequentially by using the following equations: 


where 


dx 

= g(x) 




+ 2P( t ) 

dt 




dP 

!? ? 

f p 

' ag 

dt 




ax 


. ax 

H(X) 

= 1 - 

h(x)| 


ax 


w i 


+ 2P 


ax 


h(x) } 

[H(x)Q{y<t) - h(x))]P 


. .. ( A6 ) 


. . . (A7 ) 


. .. ( A8 ) 


Here P(t) is an n x n matrix and Q is an m x m matrix. The latter matrix 
allows weighting to be assigned to the elements in the observation vector. H is an 
n x m matrix. 


By integration of equations (A6) and (A7) one can generate estimates, x(ti ), at the 
observation times, ti . An initial estimate, x(0) of the start condition, x(o), is 
required, but the estimates x(T) will usually ~ be insensitive to the choice of 
x(0), providing that T is sufficiently large. 
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TABLE 1 


Beam No. 

Nat. freq. 
(Hz) 

Damping factor 

C = y 0 */2 

Effective mass 
of journal & 
beam (kg) 

*i 

1 

33.8 

0.010 

7.73 

1.88 

2 

58.8 

0.010 

8.35 

3.53 

3 

94.3 

0.015 

9.60 

6.50 
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Figure 3. - Cross section of journal 
and beam assembly. 
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Normalized displacement Normalized displacement Normalized displacement 









Figure 5. - Averaged experimental transient 
decay results. 
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Normalized displacement 



Non-dimensional time 


(a) Variation of 6* and y* with (c) Comparison between the experi- 

non-dimensional time. mental decay curve and the theoretical 

curve, using final parameter estimates. 



(b) Comparison between the experi— (d) Comparison between the experimental 

mental decay curve and the estimated decay curve and the theoretical curve, 

displacement state variable, X(t). obtained by setting 4* = 1 and using 

the estimated value of y*. 

Figure 6. - Results of applying parametric identification to an averaged 
experimental transient decay result for € Q = 0; radial displacement. 
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Long bearing theory 


Long bearing theory 


X Shaft 1 
- Shaft 2 
o Shaft 3 


Short bearing theory 


X Shaft 1 
+ Shaft 2 
o Shaft 3 


Short bearing theory 


B rr ; radial displacement 


c rr ; radial displacement 


Long bearing theory 


Long bearing theory 


X Shaft 1 
t Shaft 2 
o Shaft 3 


\ 

Shon bearing theory 


Transverse 



Short bearing theory 

x Shaft 1 
+ Shaft 2 
° Shaft 3 

| Transverse | 


B ssi transverse displacement. 

Figure 7 . - Comparison between 
experimental results for B rr 
and B ss and the corresponding 
theoretical values, from short and 
long bearing theory. 


(*>) c ss; transverse displacement. 

Figure 8. - Comparison between exper 
imental results for m H and the cor- 
responding theoretical values from 
short and long bearing theory. 















